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Abstract 

The density crossover scaling of various thermodynamic properties of sohitions and melts of 
self-avoiding and highly flexible polymer chains without chain intersections confined to strictly 
two dimensions is investigated by means of molecular dynamics and Monte Carlo simulations 
of a standard coarse-grained bead-spring model. In the semidilute regime we confirm over an 
order of magnitude of the monomer density p the expected power-law scaling for the interaction 
energy between different chains Cint ~ P^^^^, the total pressure P ^ and the dimensionless 
compressibility = liuiq^QS{q) ~ Various elastic contributions associated to the affine 

and non-affine response to an infinitesimal strain are analyzed as functions of density and sampling 
time. We show how the size (,{p) of the semidilute blob may be determined experimentally from the 
total monomer structure factor S{q) characterizing the compressibility of the solution at a given 
wavevector q. We comment briefiy on finite persistence length effects. 

PACS numbers: 61.25.H-, 68.18.Fg, 65.20.-w 
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I. INTRODUCTION 



Compact chains of fractal perimeter. Dense polymer solutions confined to effectively 
two-dimensional (2D) thin layers are of significant technological relevance with opportuni- 
ties ranging from tribology to biology jl-3|. We focus here on the conceptionally important 
limit of self-avoiding homopolymers confined to strictly d = 2 dimensions where chain inter- 
sections are forbidden as shown by the snapshot presented in fig. [H Such systems are not 



only of theoretical 



4-|6| and computational [? 



16| but also of experimental interest 17H24|. 



n 



19 



24| . It is now 



especially since conformational properties may directly be visualized 17|, 
generally accepted J-Iisl. that at sufficiently high monomer density p and chain length 



the chains adopt compact conformations, i.e., the typical chain size R scales as 



R ^ (N/pY where u = l/d= 1/2. 



25| 



(1) 



We stress that eq. ([T]) does not imply Gaussian chain statistics since other critical exponents 
with non-Gaussian values have been shown to matter for various experimentally relevant 



)roperties 



BE, 



13| . It is thus incorrect to assume that excluded- volume effects are screened 



10l | as is approximately the case for three-dimensional melts 26| . Interestingly, compactness 
does not imply disklike shapes minimizing the chain perimeter length L. In fact the perimeter 
is found to be fractal with [gI, [121 15| 



Nn 



int 



Rdp _ N^p/d ^ith dp = d- 02 = 5/4 



(2) 



where riint stands for the fraction of monomers interacting with other chains and the fractal 
line dimension dp is set by Duplantier's contact exponent 6*2 = 3/4 B. 

Semidilute regime. Duplantier's predictions obtained using conformal invariance B i"6ly 
on the non-intersection constraint and the space-filling property of the melt. Obviously, high 



densities are experimentally difficult to realize for strictly 2D layers [17|, ll8[ since chains tend 
either to detach or to overlap, increa sing thus the number of layers as demonstrated from 
the pressure isotherms studied in ref . 18| . It is thus of some importance that eqs. ( |T|2|) have 
been argued to hold more generally for all densities assuming that the chains are sufficiently 
long 15 1 . Following de Gennes' classical density crossover scaling {4] polymer solutions may 
be viewed as space- filling melts of "blobs" of size S,{p) containing g{p) ^ p^'^{p) monomers 
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FIG. 1: We consider numerically thermodynamic properties of self-avoiding polymers in strictly two 
dimensions without chain intersections focusing on the experimentally relevant semidilute regime. 
The snapshot has been obtained for our coarse-grained model at number density p = 0.125 and 
chain length = 2048. The numbers refer to a chain index used for computational purposes. 
The large dashed circle represents a hard disk of uniform mass distribution having a radius of 
gyration Rg equal to the semidilute blob size ^ ~ 31 of the given density. While the chains adopt 
compact configurations on larger distances r ^ as characterized by the exponents v = 1/d and 
62 = 3/4, the swollen chain statistics [vq = 3/4, 02, = 19/12) remains relevant on smaller scales 
as shown in ref. [l5]. The few small open circles (not to scale) correspond to monomers interacting 
with monomers from other chains. As shown in sect. IIIIB|, the interaction energy eint due to such 
contacts remains small, even at higher densities, if the chain length is large. 



with 



e(p) ^ hg'^^ ~ l/p3/2 and 

g{p) ^ l/{htpf'^^^'-'^ ~ II p\ (3) 

Here z/q = 3/4 stands for Flory's chain size exponent for dilute swollen chains in c? = 2 
dimensions (4| and 6o = liniAr^oo -Ro(^)/^'^" for the corresponding statistical segment size. 
(Throughout the paper the dilute limit of a property is often characterized by an index 0.) 
Focusing on conformational properties it has thus been confirmed numerically that, e.g., 
nint ~ l/N""^^ = l/A^3/8 densities if > g{p) while nj^t ~ l/A^^«^2.o = for 

N g{p) with 02 = 19/12 being the contact exponent in the dilute limit. We remind that 
by matching the power laws for the dilute and dense limits at N/g ^ 1 it follows that in the 



semidilute density regime [15 1 



-.nt - ^ X (iV/,)--^ ~ p-/«iV-/«, (4) 



i.e. the fraction of monomers in interchain contact increases strongly with density^ 



Focus of present study. Having received experimental attention recently 17 



18 



20 



21 



24| . the aim of the present study is to discuss the density dependence of various thermody- 



namic properties such as the pressure P or the compression modulus K of the solution [27H29| 
focusing on the experimentally relevant semidilute regime. Comparing various numerical 
techniques we will confirm, e.g., that the dimensionless compressibility gT{p) = k^Tp/K 
scales as the blob size g{p) as one expects according to a standard density crossover scal- 



ing {4]. This is of some importance since due to the generalized Porod scattering of the 
compact chains the intrachain coherent structure factor F{q) has been shown to scale as 
NF(q) ^ A^^/(gi?)^'^~'^p in the intermediate regime 1/R ^ q ^ 1/^ of the wavevector q 



12 



15| . It is hence dangerous to determine the blob size by means of an Ornstein-Zernike 
fit as done, e.g., in ref. (171] . 

Outline. We begin the discussion in sect. Hi] by summarizing our coarse-grained model 
and the computational schemes used. Our numerical results are presented in sect. IIIII 
We remind first in sect. IIII Al the scaling of the chain and subchain size already presented 



elsewhere 



15| and characterize the size of the semidilute blob. We present then different 



energy and pressure contributions, the dimensionless compressibility g^ (sect. IIIID|) and 



the elastic moduli ?7b and rj-p characterizing, respectively, the affine linear response to an 



external homogeneous strain and the counteracting stress fluctuations (sect. IIIIE|) |30N33|. 
Some complementary information concerning the elastic Lame coeflicients A and fi [27] and 
their various contributions is referred to the Appendix. How the blob size may be determined 
in a real experiment using the scaling of the total structure factor S{q) is shown in sect. lIIIFl 
We conclude the paper in sect. IIVI where we comment on the relevance of our flndings for 
polymer blends conflned to ultrathin slits and the influence of a flnite persistence length. 



II. COARSE-GRAINED POLYMER MODEL AND COMPUTATIONAL ISSUES 



Effective Hamiltonian. The aim of the present study is to clarify universal power-law 
scaling predictions in the limit of large chain length and low wavevector q where the 
specific physics and chemistry on monomeric level is only relevant for prefactor effects 4, 
34| . As in our previous studies |l2l-ll6| we sample solutions and melts of monodisperse, 
linear and highly flexible chains using a version of the well-known Kremer-Grest (KG) bead- 

nn 

spring model [36|, |37[ . The non-bonded excluded volume interactions between the effective 
monomers are represented by a purely repulsive (truncated and shifted) Lennard- Jones (LJ) 
potential 35, Q 



u 



^b(r) = 4e [{a/rY^ - [o/rf] + e for r/a < 2^1^ 



(5) 



and u^hij) = elsewhere 



39| . At variance to the_standard KG model (i) the LJ potential 



is assumed not act between adjacent monomers '40!] and (ii) these bonded monomers are 
connected by a simple harmonic spring potential |4ll ] 

1 



-Mr - 



(6) 



with a spring constant ky^ = 676e and a bond reference length k = 0.967cr calibrated to the 
"flnitely extendible nonlinear elastic" (FENE) springs of the original KG model 36|. No 
additional stiffness term has been included and, at variance to most experimental systems 



17 



19 



23| . our chains are flexible down to monomeric scales. Effects of flnite persistence 



length are only briefly alluded to in the outlook presented at the end of the paper. 

Units and non-intersection constraint. The monomer mass m, the temperature T and 
Boltzmann's constant are all set to unity, i.e. /3 = l/k^T = 1 for the inverse temperature, 
and LJ units {e = a = m = 1) are used throughout the paper. The parameters and settings 
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FIG. 2: Normalized radial bond length distribution p{r) vs. r — for N = 2048. For all ensembles 
we obtain essentially p{r) oc exp(— /3Mb('')) as expected if the excluded volume interactions were 
rigorously switched off. Note that for MD simulations with 6t = 0.01 the tails of p{r) are slightly 
larger at low densities (solid line) compared to the MC simulations (circles) as shown for p = 0.0625 
where I = 0.9700. At larger densities the excluded volume interactions reduce the distance between 
bonded beads as shown for p = 0.875 (triangles) where I = 0.963. 

of our model make chain intersections impossible. (It has been explicitly checked that such 
a violation never occurs.) This is illustrated in the snapshot presented in fig. [T]for chains 
in the semidilute density regime at number density p = 0.125, chain length = 2048, 
chain number M = 48 and linear box dimension Lbox ~ 887. We simulate thus strictly 2D 
self-avoiding walks as required. 

Configuration sampling. Taking advantage of the public domain LAMMPS implemen- 



tation 



371] the bulk of the presented results for the more time consuming higher densities 



42] has been obtained by MD simulation integrating the classical equations of motion with 



the Velocity-Verlet algorithm 



35 



using the standard time step 6t = 0.01 



The constant temperature T = 1 is imposed by means of a Langevin thermostat [35 



4, 


16 




36] 


35. 


36, 


38 



with a friction constant 7 = 0.5. Note that the Langevin thermostat is believed to stabilize 



the integration at a larger time step (an improvement of about a factor 10 is reported) 
than necessary for microcanonical simulations [43] . We emphasize that the strong harmonic 
bonding potential — used to avoid chain intersections — corresponds to a small oscillation 
time Tb = lix ^ m I ~ 0.24. Since Xb is only an order of magnitude larger than this begs 
the question of whether configurations of correct statistical weight have been sampled. 

In order to crosscheck our results we have in addition performed Monte Carlo (MC) 
simulations which (by construction) obey detailed balance 38|, i.e. produce an ensemble 



of configurations with correct weights. A mix of local monomer moves (with displacement 



attempts uniformly distributed in a dis 
moves along the chain contours is used i 



i of radius 5r = 0.1) and global slithering snake 



34| . The latter slithering snake moves turn out to 



be efficient for exploring the configuration space at low densities up to p ^ 0.25. As one 
expects, for larger densities the acceptance rate of the snake moves deteriorates, since it 
becomes too unlikely to find enough free volume to place a new chain end. 

The comparison of ensembles generated with both methods shows that all sampled prop- 
erties are essentially identical. This can be seen, e.g., in fig. |2]for the bond length distribution 
p(r) for chains of length = 2048. Note that J dr(27rr)p(r) = 1 and J dr(27rr)p(r)r^ = 
The expected weak density effect for the bond length is clearly seen in the figure. See Ta- 
ble |T] for the root-mean square bond length /(p). Only closer inspection reveals that the MD 
method with bt = 0.01 yields slightly too large / at lower densities and the tails of p{r) are 
weakly enhanced when compared to our MC results. We shall see in sect. IIII Cl that this tiny 
numerical effect matters for the computation of the pressure in the dilute density regime. 

Parameter range. Some relevant conformational and thermodynamic properties are 
summarized in Table [T]for our main reference chain length A^ = 1024. As in ref. 15| where we 
have discussed various conformational properties of semidilute solutions and melts we scan 
over a broad range of densities p. All properties reported for p < 0.25 have been sampled 
using MC 42|. The highest density we have computed is p = 1.0 for a chain length A^ = 64. 
In cases where chain length does not matter this data set is often presented together with 
data obtained for A^ = 1024 at lower densities. Note that our largest chain length A^ = 2048 
is about an order of magnitude larger than in previous computational studies of semidilute 
solutions and melts in two dimensions 7H11| . To avoid finite system size effects the periodic 
simulation boxes contain at least M = 96 chains for the higher densities. Even more chains 
are sampled for shorter chains. 



III. NUMERICAL RESULTS 



A. Chain and subchain size 

Dilute reference limit. A density crossover scaling study requires the precise characteri- 
zation of the dilute limit Q]. Using slithering snake MC moves we have thus sampled single 
chain systems (M = 1) with chain lengths up to = 8192. To avoid the self-interaction of 
the chains with their periodic images huge simulation boxes (Lbox = 10^, p ~ 0) are used. 
As one expects j^, the typical chain size is seen to increase with a power-law exponent 
u = uq = 3/A (not shown). For the root-mean-square chain end-to-end distance Re{N) we 
obtain Rf,{N)/N'^° — )■ 6e,o ~ 0.98 for the largest chain lengths probed. The corresponding 
effective segment size for the dilute radius of gyration Rg{N) is found to become 

lim RJN)/N^'' = 6g,o ^ 0.37. (7) 

Dropping the index "g'" the latter length scale Bq is used below to make various properties 
dimensionless, allowing thus a meaningful comparison with experiments or other computa- 
tional models. 

Scaling for finite densities. That sufficiently long 2D polymer chains become indeed 
compact for all densities, as stated by eq. ([I]), is reminded in fig. [3l We present here the 
root-mean-square end-to-end distance Re{s) and the radius of gyration Rg{s) of subchains 
of length s = |m — n + ll < N = 1024 as sketched in the figure. The averages are taken 
over all pairs of monomers (n, m = n + s — 1) possible. Averaging only over subchains 
at the curvilinear chain center {n,m ^ N/2) slightly reduces chain end effects; however, 
the difference is negligible for the large chains we focus on. The limit s = N corresponds 
obviously to the total chain size which is represented in the figure for chains of length 

= 1024 (stars). As expected, the typical (sub) chain size increases with an exponent 
= z/q = 3/4 in the dilute limit (dashed lines) and with i> = 1/d for larger (sub)chains 
and densities (x ^ 1) in agreement with various numerical [7-91 and experimental studies 
[171 . 1191 ]. The subchain size is represented here to remind that not only the total chain 
becomes compact but in a self-similar manner the chain conformation on all scales 13|. 

Operational definition of blob size. In agreement with the standard density crossover 
scaling, eq. ([3]), the axes of fig. |3]have been made dimensionless by plotting y = R{s)/C,{p) 




x = (s-1)/g(p) 

FIG. 3: Root-mean-square end-to-end distance Re{s) and radius of gyration Rg{s) for subchains 
ois = m — n + l<N monomers with N = 1024. The total chain end-to-end distance Rc{N) and 
radius of gyration Rg{N) are indicated by stars (s = N). Plotting Rc{s)/S,{p) and Rg{s)/^{p) vs. 
(s — l)/g{p) yields a perfect data collapse for all densities. 

as a function oi x = {s — l)/g{p) where we define 

^ = 60^''° ~ 1.35/p^/2 and 

g = 0.09/(6[Jp)i/("»''-^) ^ 5.43/p^ (8) 

The (shghtly arbitrary) prefactors have been fitted using Rg{s) for = 1024 and = 2048 
for the densities p = 0.0625 and p = 0.125. This yields a perfect data collapse especially 
considering that a broad range of densities is considered. Please note that the asymptotic 
power-law slopes for x <^ 1 (dashed lines) and x ^ 1 (bold lines) for the radius of gyration 
intersect exactly at {x,y) = (1, 1). Since Re{s) becomes compact more rapidly than -Rg(s), 
a blob size defined using Re{s) would be slightly smaller. 

Scaling for large densities. Interestingly, the blob scaling of the (sub)chain size is even 
successful for our highest densities where the blob picture clearly breaks down for some other 
properties as discussed, e.g., in sect. IIII CI below. This is due to the fact that R{s) is set in 
this limit by the typical distance dcm ~ i-^/pY^'^ between the chain or subchain centers of 




FIG. 4: Various mean energy contributions per monomer vs. density p for chain length N = 1024: 
Total potential energy etot, bonding energy Cb, excluded volume energy Cnb and interaction energy 
Cint between monomers of different chains. The power-law slope predicted for the semidilute regime 
is indicated by the bold line. The vertical arrow indicates p* for N = 1024 as determined from the 
pressure isotherm, eq. p3|) . discussed below. Non-universal behavior due to the specific monomer 
interactions is seen at higher densities with p > p** ~ 0.5. 

mass and this irrespective of the physics (monomer size, persistence length, . . .) on small 
scales. Since p ^ g/C,'^ is imposed, a different choice of the blob size only leads to a shift of 
the data along the bold power-law slope u = 1/d. Only for sufficiently low densities where 
(sub) chains smaller than the blob can be probed it is possible to test the blob scaling and 
to adjust the prefactors. 

B. Energy 

From the numerical point of view the simplest thermodynamic property to be investigated 
here is the total mean potential energy per monomer etot due to the Hamiltonian described 
in sect. HIl As shown in fig. HJ it is essentially density independent and always dominated 
by the bonding potential e^- Due to the harmonic springs used we have ~ k-BT/2 as 




FIG. 5: Interchain interaction energy per monomer eint vs. chain length N for different p. The 
dashed and sohd power-law slopes with — 1^02 = —19/16 and —3/8 represent, respectively, the 
expected asymptotic behavior for dilute and dense solutions. 



expected according to the equipartition theorem [4l|. Noticeable (albeit weak) corrections 
to this value are only found for our highest densities where bonded monomer pairs are pushed 
together by the excluded volume interactions as already seen in fig. [2l The total non-bonded 
excluded volume interaction per monomer e^h becomes constant at low densities where it is 
dominated by the excluded volume interaction of curvilinear neighbors on the same chain. 
The non-bonded energy Cnb increases for larger densities, but remains always smaller than 
the bonded energy e^. Values of Cnb for = 1024 are included in Table [B 

From the theoretical point of view more interesting is the contribution to the total ex- 
cluded volume interaction due to the contact of monomers from different chains measured 
by Cint- For not too high densities Cint is expected to scale as the fraction nint of monomers 
in interchain contact mentioned in the Introduction. As indicated by the dashed line the 
interchain energy eint is proportional to the density in the dilute regime (dashed line) due 
to the mean- field probability that two chains are in contact. At higher semidilute densities 
above the crossover density p*{N) up to p** ~ 0.5A^° a much stronger power-law exponent 
~ 21/8 is seen in agreement with the established scaling for njnt, eq. (jlj). The observed 



p-dependence of Cmt in the semidilute regime is thus traced back to the known values of 
the universal exponents v and Q2 in the dilute and dense limits. The interaction energy 
Cint is seen to increase even more strongly for densities p > p** where the semidilute blob 
picture becomes inaccurate. At variance to the other mean energies, eint has a strong chain 
length effect as revealed in fig. [51 The indicated power-law slopes correspond to the expected 
exponents z/6'2 = 19/16 and 3/8 for, respectively, the dilute (dashed line) and dense (bold 
lines) density limits. Since the A^-scaling does not require the existence of sufficient large 
semidilute blobs, it even holds for our highest "melt" densities. 



C. Pressure 



Definitions and virial equation. While the energy contributions discussed above cannot 
be probed in a real experiment, the osmotic pressure of 2D polymer systems can be accessed 



experimentally for polymer solutions at the air- water interface 18|, l20l-l24l]. As usual for 
pairwise additive interactions the mean pressure P = Pid + Pcx is obtained in our simulations 
as the sum of the ideal kinetic contribution P^^ = k^Tp and the excess pressure contribution 

Pe. = Pnb + Pb = (W) /V. (9) 

Here, Pnb stands for the non-bonded excess pressure contribution, Pb ioT the bonded pressure 
contribution, V for the d-dimensional volume, i.e. the surface L^^^ of our periodic simulation 



box, and W for the internal virial 



35 



,l38| 



I I 

with fi being the force of the interaction I between two beads i and j at a distance ri = 
I IZLj "ZLi I I and w{ri) = riu'{ri) the virial function associated with the bonded and non-bonded 



pair potential u{ri). (The sum Y2i stands for the double sum Yli<j '^^^^ monomers 
of the solution.) 

Pressure contributions. We present in fig. [6] the different contributions to the total pres- 
sure P = Pid -\- Pnb ~l~ Pb clS a function of density p focusing on the chain length = 1024. 
The vertical axis is rescaled by a factor /3/p to present the contributions per monomer as 
in fig. m for the different energy contributions. The dilute limit, where PP/p — )■ 1/A^, is 
indicated by the thin dashed line, the semidilute limit by the bold line representing the 
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FIG. 6: Various contributions to the isothermal pressure for chains of length N = 1024 as a 
function of density p: the total pressure P = Pid + Pex, the ideal kinetic pressure = k^Tp, the 
(negative) bonded pressure contribution —Pb, the non-bonded LJ interaction contribution Pnb, the 
(negative) excess pressure —Pex = —{Ph + -fnb) and the pressure contribution Pint due to the LJ 
interactions between monomers on different chains. The thin dashed line indicates the dilute total 
pressure limit P/3/p = 1/N, the bold line the exponent duQ/{duQ — 1) — 1 = 2 according to eq. (jlip 
and the dashed bold line the exponent expected for Pint in the semidilute regime. 

expected power law PP/p ~ discussed below. The rescaled non-bonded pressure Pnb 
becomes constant for densities p < p** ~ 0.5 with a plateau value PnhP/p ~ 0.23 due to the 
intrachain interactions of closely connected neighbors along the chain as the corresponding 
energy contribution Cnb presented in fig. HI 

In the same density regime we have PhP/ p ~ —1.23 for the pressure contribution Pb due 
to the bonding potential. We remind that for asymptotically long non-interacting phantom 
chains one would have Pb/3/p = —1 irrespective of the details of the bonding potential 
44l |. At variance to the energy contribution Cb the corresponding pressure contribution Pb 
is thus slightly changed by the presence of the excluded volume potential. The pressure 
contribution Pb must be more negative, of course, since the excess pressure Pex = -Pnb + -Pb 
must essentially cancel (for very large chains) the ideal pressure Pid indicated by the thin 



solid line. Note also that for our highest densities where the bonded monomers are pressed 
together both Pb and Pex become positive, i.e. are not represented in the figure. 

The pressure contribution PmtP/ p due to interchain monomer contacts (squares) shows, 
not surprisingly, the same power-law exponents as in the interchain interaction energy Cint 
discussed above. In the semidilute regime we confirm P\nt(^/ p ~ p^^^^ as indicated by the 
bold dashed line. The expected dilute limit P-mtP/ p ~ p (dash-dotted line) is unfortunately 
not yet confirmed numerically due to insufficient statistics and the error bars (not given) 
become much larger than the symbol size {45 1. 

Scaling with chain length. The total pressure P is presented in fig. [7] for a broad range 
of chain lengths A^. Obviously, P increases strongly with density p and the chain length 
only matters in the dilute limit for short chains or small densities where ultimately the 
translational entropy of the chains governs the free energy j4| as indicated by the thin dashed 
horizontal lines. This (theoretically trivial) limit turns out to be numerically challenging 
since a large positive term, the ideal pressure contribution Pjd, is nearly canceled by a large 
negative term, the excess pressure Pex ~ Pb, which requires increasingly good statistics as 

becomes larger. The precise determination of P becomes surprisingly time consuming 
even if the configurations are sampled by means of slithering snake MC moves (sect, ini). 
Given sufficient numerical precision this yields, however, the expected dilute pressure as is 
seen from the main panel of the figure. This is different if systems are computed using MD 
with the standard time step 5t = 0.01 as shown in the inset of fig. [71 The bond oscillation 
time Tb/lOO is indicated by the vertical arrow. As we have already seen in fig. [21 the bonds 
become slightly stretched if 6t ^ Xb/lOO, i.e. they become too tensile which corresponds to 
a too negative Pb. Unfortunately, the bonding potential is that strong, i.e. Tb/lOO so small, 
that it gets too time consuming to compute the correct pressure in this density limit using 
MD [46]. 

Semidilute density regime. Of experimental interest beyond these computational issues 
is the intermediate semidilute density regime indicated by the bold power-law slope. The 
universal exponent can be understood by an elegant crossover scaling argument given by de 
Gennes where the pressure is written as /3P = p/N x f{p/p*) with p*(N) ^ N/R'^{N) 
crossover density and f{x) a universal function 



47j . Assuming P to be 
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FIG. 7: Total pressure P/3/p as a function of density p for different chain lengths A^. Unconnected 
LJ beads (N = 1) are indicated by the small filled circles. The data for iV = 16 is taken from the 
instantaneous pressure obtained by MD simulation increasing extremely slowly the box size. The 
computation of larger chain lengths is much more delicate. The thin dashed lines give the osmotic 
pressure in the dilute limit for = 1, 16, 64, 256 and 1024 (from top to bottom). The pressure 
becomes A^-independent with increasing N and p. The statistics deteriorates for small p and large 
N. The bold line indicates the power-law exponent dvQ/{dvQ — 1) = 3 expected in the semidilute 
regime according to eq. (Ilip . Inset: Pressure as a function of the MD time step 5t for N = 1024 
and p = 0.0625. The horizontal line corresponds to our MC result. As indicated by the vertical 
arrow the correct pressure in the dilute limit is obtained only if 5t <C rt/lOO. 



chain length independent for a; 3> 1 this implies f{x) ~ x^^^'^'^° and, hence. 



(11) 



where we have used eq. ([3]) to restate the well-known relation between pressure and blob size 
^ j4|. The predicted exponent fits the data over about a decade in density where the blob 
size is sufficiently large. Additional physics becomes relevant for densities around p** where 
the LJ excluded volume starts to dominate all interactions, i.e. the pressure of polymer 
chains approaches the pressure of unbonded LJ beads (ffiled circles). 



Universal amplitude. The choice of the axes of the main panel of fig. [7] may make the 
comparison to real experiments difficult. Choosing as a (natural but arbitrary) length scale 
the effective segment size 60 associated to the dilute radius of gyration, eq. ([7]), one may 
instead plot the rescaled pressure y = pPh^ as a function of the reduced density x = b^p. 
In the semidilute regime this corresponds to a power-law slope y = apx^ with a power- 
law amplitude ap ^ 83.4. The dimensionless amplitude ap (or similar related values due 
to different choices of Bq) should be compared to real experiments or other computational 
models. As long as the blob size is sufficiently large, i.e. p small enough, molecular details 
should not alter this universal amplitude. Persistence length effects, e.g., change 60 but not 
ap. Using eq. this implies 

(3P ^ 0.16ap/^2 ^ 13.7/^2 ^^2) 

with universal prefactors (within the operational definition based on the radius of gyra- 
tion) allowing thus to determine the blob size from the experimentally obtained pressure 
isotherms. We note finally that by matching the dilute asymptote with the semidilute pres- 
sure regime the prefactor of the crossover density p* may be operationally defined as 

p* = a^-"^^" X N/R'^^ ^ 0.845/A^^/2 (13) 

with Rg being the radius of gyration in the dilute limit. This implies p* ~ 0.02 for our largest 
chains with = 2048. Considering that the semidilute regime breaks down at p ^ p** this 
limits the semidilute scaling to about an order of magnitude in density. 



D. Compressibility 



Definitions. Being an isotropic liquid, a polymer solution is described in the hydro- 



dynamic 
modulus 



imit by only one experimentally relevant elastic modulus, the bulk compression 



28 



29| 



dP 

K =1/kt = ksTp/gT^^ = p—, 



(14) 



with kt being the standard isothermal compressibility and (^t.n the "dimensionless com- 
pressibility" for systems of finite chain length A^. We use here the additional index A^ 
to distinguish (yfT,N from the dimensionless compressibility for asymptotically long chains 
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FIG. 8: Determination of dimensionless compressibility 5t(p) = liniTv-^oo 5t,n(/o): (a) (7t,n(/o) 
obtained by various means for several chain lengths N. The thin lines correspond to a polynomial 
fit to P{p) for = 1 (lower thin line) and the largest chains available (upper thin line). The large 
circles have been obtained using the Rowlinson formula, eq. (|16p . for N = 1024. All other data 
correspond to the plateau of the total structure factor S{q, N) for small wavevectors q. The bold 
line indicates eq. ([3]), the two horizontal dashed lines the dilute limit for = 16 and N = 256. (b) 
Data collapse of 1/(7t,n — 1/A^ for different N with 5(t,n obtained from the total structure factor. 



grp = liniAr^oo fl'T.N- Due to the translational entropy of the chains (van't Hoff's law) both 



quantities are related by 



26| 



1 



fl'T,N 



1 

N 



1 



(15) 



i-e. Qt^^ and are expected to differ strongly for small densities where must be large. 
Our aim is to determine (7t,n precisely comparing different techniques and to extrapolate 
then using eq. (IT5]) to g^; which should scale in the semidilute regime as the number of 
monomers per blob, eq. (|3]). 

Measurements. Using eq. ( !T^ the bulk modulus and/or the dimensionless compress- 
ibility can be obtained, of course, from the pressure isotherms discussed in the previous 



subsection. This is best done by fitting a spline to y = log(/3P) as a function of x = log(p). 
Tlie resulting curves for = 1 and the largest chain lengths available are represented by 
the thin lines in the main panel of fig. E] where gT,^ is traced as a function of density p. A 
disadvantage of this method is of course that P{p, N) must be known for a large number of 
densities, especially for large p where the pressure increases strongly. 

Alternatively, the compressibility for one specific density can be directly computed using 



the Rowlinson stress fluctuation formula (28 



35| 



K = P + 7]B-VF (16) 

where we emphasize that K depends explicitly on the total pressure P. The second contri- 
bution r^B represents the so-called "hypervirial" 

where / stands again for a pair of monomers i < j. The index B indicates that this contribu- 
tion corresponds to the famous Born approximation of the elastic moduli of solids assuming 



affine displacements under an applied infinitesimal homogeneous strain [30H32| . (The used 
notation becomes transparent from the exact relation, eq. (lA7p . given in the Appendix.) This 
affine approximation overpredicts the free energy change in general. The overprediction of 
the compression modulus is "corrected" by the excess pressure fluctuation 

77F = f3V {{{W/Vf) - {W/Vf) > 0. (18) 

Results obtained using eq. f|T6l) and local MC moves are indicated by the large circles in 
fig. [HI While the computation of K using the Rowlinson formula is straightforward for 
dense systems (p > 0.5), this becomes for numerical reasons more and more difficult with 



decreasing density, as further investigated in sect. IIIIE _48|. 



The bulk of the data presented in fig. [8] stems from 



gT,^ = p/PK = \imS{q,N) (19) 

(J— 5>0 

from the plateau in the low-wavevector limit of the total structure factor S{q, N) which is 
further discussed in sect. lIIIFl (Equation (fT9|) takes advantage of the fact that our monomers 
are indistinguishable. It would become much more intricate if the beads were, e.g., polydis- 
perse j29|-) Due to our large box sizes (especially in the low-p limit) this method provides 
over the whole density range reliable numerical values only requiring the analysis of about 
1000 more or less independent configurations. 



Interpretation of data. In agreement with eq. f|T4|) and the scahng of the pressure dis- 
cussed above, the dimensionless compressibihty (7t,n is seen in fig. [8] to be strongly A^- 
dependent for short chains and small densities where the translational entropy matters. 
This A^-dependence is, however, perfectly characterized by eq. (JTSil as explicitly verified by 
the scaling collapse presented in the inset. The scaling being successful for even rather 
small chains, this allows us to determine even using short chains the large- limit g^^ for all 
densities (indicated by stars). Consistent with eq. (fTTj) the bold line indicates the power-law 
asymptote 

for the semidilute regime. Using eq. ([H]) this implies g ~ 23(7t for the number of monomers 
per blob which may thus be obtained from the dimensionless compresssibility. 

E. Stress fluctuations 

Introduction. Using different numerical techniques we have determined in the preceding 
subsection the compression modulus Kf3/p = 1/(?t,n relevant for real experiments. As above 
for the energy and the pressure we discuss now in more detail various numerically accessible 
contributions to K. One aim is to better analyze the already mentioned numerical difficulties 
encountered when probing the compressibility for densities below p** using eq. ( TTGl) . Our 
second aim is to characterize the fiuctuations of the (excess) pressure considered in sect. lIII Cl 
Various elastic moduli are shown in fig. |9] as functions of the density p and in fig. [10] and 
fig. [11] as functions of the sampling time t. All data presented here have been obtained by 
MC simulation, i.e. there is no discretization problem as for MD, and the time is given 
in MC Steps (MCS) of the local monomer jump attempts Issl. The vertical axes of the 
figures are made dimensionless by a factor /3/p as in fig. H] and fig. [6] i.e. we focus on 
the elastic free energy contribution per monomer. While we have insisted above on the 
A^-dependence due to the translational entropy of the chains, eq. f[T5]) . we consider now such 
high densities and/or large chains that this (additional) complication may safely be ignored 
for the presented data. 

Density dependence of contributions to K . The compression modulus K{p) obtained 
using the Rowlinson formula for chains of length A^ = 1024 replotted in fig. [9] may indeed 
be considered as A^- independent. Also given are the (rescaled) contributions to K according 
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FIG. 9: Contributions to K for N = 1024 vs. density p: Compression modulus K = P + rj^ — rjp, 
pressure P, hypervirial r/e and fluctuation contribution rjp. Also given are the hypervirials r/B,nb 
and 7yB,b and the pressure fluctuation contributions r/F,nb and ryF,b associated to the non-bonded and 
the bonded interactions, respectively. Filled symbols refer to the hypervirial (Born) contributions. 
All data points for p = 1.0 have been obtained for = 64. The horizontal line corresponds to 
^B,b = /3^b^^/4 as expected for non-interacting chains, eq. (|22p . and the bold solid line to the 
compression modulus in the semidilute regime obtained using eq. (j20p . The contribution ?yF,seif,nb 
associated to the two-point or self correlations of the non-bonded interactions is represented by the 
dashed line. 

to eq. (fTB]) : the mean pressure P, the hypervirial rj^ and the excess pressure fluctuation 
rjp. While the pressure is small over the full density range, r]Bf3/p and rjpP/ p are seen to be 
essentially density independent and of same magnitude. A reasonable numerical estimation 
of the compression modulus K thus requires a precise determination of its two leading 
contributions. Note that the logarithmic representation masks the noise naturally present 
in the data, especially for lower densities. (This may be better seen from the related data 
presented in the inset of fig. [15] given in the Appendix.) Since with decreasing density K 
becomes rapidly orders of magnitudes smaller than r/B ~ ?7f; it is not surprising that we 
have not been able to obtain reliable values for K from eq. (fT6|l below p ~ 0.25. 



Moduli associated to different interactions. Applying eq. f lT7|) and eq. f|T8l) to the non- 
bonded and bonded potential contributions one obtains the hypervirials ?7B,nb and T^B.b and 
the stress fluctuations r^p.nb and ?7F,b- Since the hypervirial is linear with respect to the 
different interactions we have r/B = VB,nh + VB,h, while rjp = ?7F,nb + VF,h + Vfm^ where the 
last term r^p.mix characterizes the correlations between the bonded and non-bonded stresses. 
(Being always much smaller than the other contributions r/p.mix is not discussed here.) As 
can be seen from the figure we have for all densities 

VB^VF^ VB,h ~ r]F,h > ?7B,nb ~ %,nb, (21) 

i.e. the bonded contributions to ?7b and r]p dominate numerically by far the non-bonded 
interactions. (Obviously, this does not mean that the difference r^B.nb — VF,nh is irrelevant 
compared to the difference ?7B,b — VF,h-) The plateau value indicated by the thin horizontal 
line is thus readily computed from the hypervirial 

r/B,b/3/p = ^ (r(2r - h)) ^ (22) 

with r being the distance between bonded monomers. In the second step of eq. ( l22l) it was 
used that = (r^) ^ (r) /b ~ ll due to the stiffness of the bonding potential. 

Two-point stress correlations. Why is a similar value expected for the fluctuation con- 
tribution rjp ^ i]F,h'^ To see this we remind first that the fluctuation term, eq. f|T8|) . may be 
rewritten quite generally as 

= ^E(H^')'>-H^')n (23) 

i.e. the fluctuation contribution to K contains not only two-point correlations (/ = /') but 
also three- and four-point correlations (/ ^ I'). If we assume that to leading order only the 
two-point or "self correlations matter for the bonded interactions it follows that 

r/F,b/?/p ^ ^ {{Pwirf) - {Pwir)f) (24) 



d 



\{r-k)V) (25) 



(26) 




t [MCS] 

FIG. 10: Elastic properties for p = 0.75 and N = 1024 as functions of the sampling time t: 
Compression modulus K, total pressure P, hypervirial rye, total excess pressure fluctuation t/f and 
the contributions r/F,nb and T^F.b due to the non-bonded and bonded interactions. While P and r/B 
are (essentially) constant, r]p{t), r/F,nb(i) and ?/F,b(i) are seen to converge only slowly with sampling 
time and, hence, the compression modulus K{t). 

where r stands again for the length of a bond. In the second step the small squared excess 
pressure term {f3w{r))'^ is neglected and we have finally taken (again) advantage of the 
stiffness of the harmonic potential. The self- or two-point correlation contribution to rj-p 
can of course be computed directly. For the total contribution r^p.seif and the contribution 
'7F,seif,b of the bonded interactions we obtain essentially eq. ( I26l) for all densities. This is not 
represented in fig. [9] since these values could not be distinguished (in the double logarithmic 
representation chosen) from the values rj^, rjp already given. Instead we show the self- 
contribution r^F.scif.nb associated to the non-bonded potential (dashed line) which is seen to 
be essentially identical to ?7B,nb ~ VF,nh [49] . 

Time dependence of the compression modulus. We have seen above that one difficulty 
to determine the compression modulus for dilute and semidilute polymer solutions using 
the Rowlinson stress fluctuation formula stems from the fact that a large hypervirial rj-B is 
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FIG. 11: Monotonous decay of the compression modulus K{t) as a function of the sampHng time 
t for different densities p for one chain length N = 1024. Inset: The smaher the density, the 
later K{t)f3/p levels off as indicated by the solid horizontal lines. Main panel: Scaling plot of 
K{t)/K{t ^ oo) as a function of the reduced MSD {h{tY/'^ / i{p)Y with i{p) being the blob size 
assuming eq. dH). The scaling is not perfect, especially not for larger densities, but catches the 
essential density effect. 

essentially compensated by an equally large stress fluctuation term rj^. This requires a high 
precision for determining both contributions. The computational difficulty is in fact not 
Tj-Q which is readily obtained to high precision, but the fluctuation contribution rj-p which is 
found to require a substantial subvolume of the configuration space to be sampled. This 
point is corroborated in fig. [10] presenting various (ultimately static) elastic contributions 
as functions of the sampling time t for density p = 0.75. The presented data have been 
obtained exclusively using MC simulations with local monomer displacements. Using time 
series where instantaneous properties relevant for the moments are written down every 10 
MCS. All reported properties have been averaged using standard gliding averages jssl, i.e. 
we compute mean values and fiuctuations for a given time interval [to, ti = to+t] and average 
over all possible intervals of length t. It is seen that simple means such as the pressure P and 



the hypervirial r/B reach immediately their asymptotic values. The open circles correspond to 
K computed according to eq. ( fT6l) for a canonical ensemble at constant volume. Measuring 
fluctuations rj-p must vanish if only one configuration is measured (t = 0). Since riY{t) 
increases for small times, K approaches the asymptotic value from above requiring about 
10^ MCS to reach the plateau. (That the fluctuation contribution rj-p^nhif) of the non-bonded 
interactions becomes similar for large times is by accident for the density presented in the 
figure.) 

Density scaling of time dependence. The time dependence of K{t) is further investigated 
in fig. [H] where we compare different densities p for = 1024. As shown in the inset K{t)(5/ p 
decreases monotonously both with decreasing density and increasing sampling time. (We 
have included here in addition a data set for p = 1.0 and = 64.) The plateau values 
K{t — )■ oo) for long times are marked by the horizontal lines. It is seen that increasingly 



more time is needed to reach the plateau if the density is lowered |48 |. 

The p-scaling of the crossover time is clarified in the main panel where we have plotted 
the rescaled vertical axis K{t)/K{t — )■ oo) as a function of the measured mean-square 
displacement (MSD) h{t) for each density. Hence, instead of the sampling time we use 
the volume h'^/'^it) explored by the beads as coordinate. The horizontal axis is then made 
dimensionless by reducing the MSD using the blob size ^{p) given by eq. (|8]). As can be 
seen we observe a satisfactory scaling collapse for all data sets. We emphasize that the time- 
independent thermodynamic limit is reached when the beads have explored the semidilute 
blob, i.e. h{t) ~ ^"^{p)- That the collapse is not perfect for our highest densities p > 0.5 is 
expected since the semidilute blob scaling breaks down in this limit. Data for lower densities 
are warranted in the future to demonstrate the suggested blob scaling unambiguously. 

F. Total monomer structure factor S{q) 

The compression modulus K discussed above describes the linear response of the systems 
in the hydrodynamic limit for wavevectors g — ?■ 0. We turn now to an experimentally 
highly relevant reciprocal space characterization of the polymer solution characterizing the 
compressibility at a given wavevector q: the "total monomer structure factor" 0, 

Siq) = - Yl <^^P h ■ (^n -rj]). (27) 

n,m=l 
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FIG. 12: Total monomer structure factor S{q) vs. wavevector q for N = 1024 and several densities 
p. In the dilute limit S{q) corresponds to the single chain form factor FQ{q) indicated by the small 
filled circles. The plateau values in the low-g limit for finite p have been used for the determination 
of 5t,n(p) as shown for p = 0.125 and p = 0.5 (thin horizonal lines). At smaller wavevectors S{q) 
is found to decay monotonously only for densities p < 0.5. At higher densities S{q) increases with 
q in qualitative difference to the RPA formula, eq. (j29p . which is indicated by the thin line for 
p = 0.125 using the directly measured intramolecular structure factor F{q). 

In our simulations all rimon monomers of the simulation box are assumed to be labeled and the 
average (...) is performed over all configurations of the ensemble and all possible wavevectors 
of length q = \q\. As already discussed in sect. 1111 Dt the isothermal compressibility of the 
solution may be obtained from the plateau of the structure factor in the low-g limit according 
to eq. ( IT9l) . The smallest possible wavevector is of course 27r/Lbox since the wavevector q 
must be commensurate with the simulation box. As can be seen from fig. [12] our box sizes 
allow for a precise determination of (^t.nIp) for all densities p > 0.0625. Only chains of 
length = 1024 are presented for clarity. Since above g > 3 the monomer structure 
becomes important for all densities and, being interested in universal physics, we focus 
below on smaller wavevectors. For comparison we have also included the intramolecular 



structure factor F{q) = -Fo(g) obtained for one single chain in a large box (p = 0). As 
emphasized by the dashed line FqIq) is characterized in the intermediate wavevector regime 
by the power-law decay 

{boqp'-^ Kg a 
with Co ~ 0.5 being a dimensionless amplitude and a representing the monomer scale. 

Returning to finite monomer densities we remind the "random phase approximation" 
(RPA) [41 

1 1 1 

- + (29) 



Siq) 9T F{q) 

with F{q) being the intramolecular structure factor at the given density. The RPA is sup- 
posed to relate — at least for not too high densities — the total structure factor S{q) to 
the dilute chain form factor F{q) ^ F^lq) for q ^ 4, |26|]. Note that gic stands for 
the excess contribution to the dimensionless compressibility in agreement to eq. (fT5|) and 
to 1/F{q) — 7- in the low-g limit. By construction S{q) is fitted by the RPA in the 

low-g and large-g limits as shown by the thin solid line for p = 0.125. Interestingly, even 
for low densities the crossover regime between both g-limits at g 27r/^ is, however, only 
inaccurately described jsO]. The monotonous decay of S{q) implicit to eq. (l29ll is indeed 
observed for the semidilute densities with p < 0.5. For higher densities, where the blob 
picture breaks down due to monomer physics, S{q) becomes essentially constant for all q up 
to the monomeric scale. 

The scaling of S{q) in the semidilute regime is further investigated in fig. [T3] where 
y = 5'(g)/5'T,N is plotted as a function of x = g,^/27r with ^ being set by the matching of the 
radius of gyration in the dilute and semidilute limits, eq. ([8]). As indicated by the vertical 
dashed line the asymptotic slopes for small and large wavevectors are found to intercept at 
X = 1 \ Using eq. ( 128|) for the dilute chain form factor this implies ^ = 6o2vr((7T/co)'^''- Since 
on the other side ^ = bog'^°, it follows that the number g of monomers spanning the blob 
is given hj g ^ 23(7t as already stated above, eq. (l20l) . Our operational prefactor setting, 
eq. dH]), thus corresponds to an experimentally (in principle) measurable choice. It suggests 
that future experimental work may proceed in a similar manner by fixing ^ (or equivalently 
g) from the matching point of both g-limits of the total structure factor rather than by 
imposing an inappropriate Ornstein-Zernike fit to the intramolecular form factor F{q) or 
the total structure factor S{q) [17]. 




X = q ^(p)/27r 

FIG. 13: Rescaled structure factor y = S{q)/gT^^{p) as a function of reduced wavevector x = 
q^{p)/2Ti. Assuming eq. ([8]) for the blob size ^, the asymptotic small-g (horizontal line) and large-g 
(bold dashed line) limits match at x ~ 1 as indicated by the vertical dashed line. 

IV. CONCLUSION 



Summary. In this paper we investigated numerically the crossover scaling of various 
thermodynamic properties of solutions and melts of a generic bead-spring model of self- 
avoiding and highly flexible polymer chains without chain intersections confined to strictly 



two dimensions. As expected js, 6, 15 1 the typical interaction energy Cint between monomers 
from different chains was shown (fig. [5]) to scale as Cint ~ with z/ = 3/4 and 62 = 19/12 

for dilute solutions and u = 1/d = 1/2 and 62 = 3/4 for sufficiently large chains with 
^ g{p) ~ 1/p^ where the chains adopt compact configurations of fractal perimeter 
dimension. In the semidilute regime we confirmed over an order of magnitude in density p the 
power-law scaling for the interaction energy Cjnt ~ p^^^^ (fig. Hj), the pressure P (fig. [7]) 
and the dimensionless isothermal compressibility (jfT ~ 1/p^ (fig. E]) expected theoretically in 
the limit of asymptotically long chains. Polymer specific technical difficulties associated with 
the numerical determination of pressure (inset of fig. [7]) and elastic moduli (figs. [QlfTT] and [T5|) 
at small and intermediate densities have been discussed. The elastic contributions r^B and 



rjp associated, respectively, to the affine and non-afiine response to an imposed homogeneous 
strain were analyzed as functions of density and sampling time. We have emphasized that 

r]B^r]F^ r]B,h ^ VF,h ~ hll/ip :^ K :^ P (30) 



for all but the highest densities 49|]. (Similar relations exist for the Lame coefficients are 
discussed in Appendix |Al) The stress fluctuations rip{t) become only time independent if 
distances corresponding to the blob size ^{p) are probed (fig. [TTll . Returning finally to 
experimentally relevant properties we showed how the size of the semidilute blob C,{p) niay 
be determined in a real experiment from the pressure isotherms, eq. f|T2|) . or the total 
monomer structure factor S{q) characterizing the compressibility of the solution at a given 
wavevector q (fig. [131) • 

Polymer blends. As already argued elsewhere 0, [isl, the presented numerical results 
for the interchain interaction energy Cint for homopolymer systems implies an enhanced 
compatibility for polymer blends confined to ultrathin films in qualitative agreement with 



recent experimental studies 2J]. We remind that the critical temperature of unmixing 
Tc should be proportional to the typical interaction energy Nci^t between different chains. 
According to eq. (jlj) and using eq. ([3]) it follows thus that Tc ~ p^i/Sjys/s compact 
chain limit. This scaling is in fact consistent with recent MC simulations of symmetrical 
)olymer mixtures using a version of the bond-fluctuation model in strictly two dimensions 
lol . The predicted exponent 5/8 = 0.625 for the A^-dependence agrees well with the 



value 0.65 obtained by fltting for all computed chain lengths [10|. The corresponding 
strong power-law increase Tc ~ p^^^^ with density has to our knowledge not been probed yet 
which calls for an experimental and/or numerical veriflcation focusing on semidilute polymer 
blends. 

Persistence length effects. Returning to homopolymer solutions we note flnally that, at 



variance to most experimental systems 



17 



19 



23| , we have assumed here that the chains are 



flexible down to monomeric scales. As long as the persistence length remains smaller than 
the blob size, this should in fact not alter the suggested scahng properties. In order to bring 
our computational approach closer to experiment we are currently computing systems with 
flnite persistence length. In addition to the coarse-grained model Hamiltonian presented in 
sect. [TTl a local stiffness potential kg{l — cos(6')) is applied with 6 being the angle between 
adjacent bonds of a chain. As shown by the pressure isotherms presented in flg. [TUfor chains 
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FIG. 14: Reduced pressure P/3/p vs. density p ioi N = 256 and several stiffness energy penalties 
kg. In agreement with previous computational findings jsi] the pressure P{kQ) increases with ke 
for low densities but decreases for high densities. Inset: Snapshot of a subvolume of a larger box 
showing local order and hairpin-like defects obtained for N = 256, p = 0.5, Lbox ^ 313 and kg = 16 
starting from an isotropic and flexible system with /cg = 0. Due to the hairpins the equilibration 
dynamics becomes sluggish and strong hysteresis effects are observed if systems with different 
dynamical pathways are compared. 



of length = 256 the scaling with respect to density for sufficiently long chains remains 
essentially unchanged. In fact, since the blob size decreases very strongly with persistence 
length, i.e. N/g increases, a finite rigidity even speeds up the convergence to the predicted 
asymptotic behavior. As expected from the decreasing blob size ^, the pressure /3P ~ l/^'^ is 
seen to increase with kg for small densities. Obviously, if the stiffness penalty kg and/or the 
density p become too large, a nematic chain alignment becomes relevant (at least locally) as 
shown by the snapshot shown in fig. [TH Due to this additional effect P is found to decrease 
with kg for large p in agreement with the numerical study by Dijkstra and Frenkel 51|. 
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FIG. 15: Lame coefficients and related properties for N = 1024. The vertical axes is made dimen- 
sionless by means of a factor (3 /p. Inset: Born Lame coefficients Ab = and stress fluctuations 
— Ap and — //p as functions of p. The data for dilute and semidilute densities is weU approximated 
by eq. ()A12p . The vertical axis being linear, the noise appears more strongly than for tjb and r]p in 
fig. [9] where a logarithmic representation is chosen. (The data point for p = 1.0 has been obtained 
for N = 64.) Main panel: Time dependent shear modulus G{t) = /^b + f^F{t) — Pcx for several 
densities as a function of MSD h(t). The data for different densities scale if a density independent 
local scale a is assumed for the horizontal axis. 
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Appendix A: Lame coefficients and related properties 



Definitions. For consistency with previous numerical work 3l|, l32|] we present here the 
different contributions to the Lame coefficients A and fi of the solution As may be seen 
by thermodynamic and symmetry considerations, the compression modulus K and the shear 
modulus G of any isotropic and homogeneous system in d dimensions may be rewritten as 



K = {X + P) + -G, 
G = jj-P 



(Al) 
(A2) 



where we follow the notation of ref. 32[. We do this to emphasize the explicit pressure 
dependence which is often (incorrectly) omitted j^^] as clearly pointed out by Birch 52 1 and 
Wallace [ssl- Note that K = X + fimd = 2 dimensions. Since by definition of a liquid the 
shear modulus must vanish for our systems, G = 0, the Lame coefficient /z is simply given 
by the total pressure P (as indicated in fig. [9]). Hence, 



K = Pid + Pcx + Ab - (-Af) 

= G = /iB - (-/ip) - ^ex 



(A3) 
(A4) 



where we have rewritten the Lame coefficients as 



32| 



A = Ab — (— Af), 
fi- Pid = fJ'B- (-/^f)- 



(A5) 



Please note that the only contribution due to the kinetic energy of the particles is contained 
by the ideal gas pressure Pjd indicated for /i. Kinetic energy contributions to the elastic 
moduli are removed as far as possible in view of the fact that MC results are considered 
here. 

Born Lame coefficients. The ffist contributions indicated on the right hand-side of 
eq. ( 1A5I1 are the so-called "Born Lame coefficients" 



riu'{ri)) 



(A6) 



where the index / stands for the interaction between two monomers i < j, ri for the length 
of the vector between both monomers and xi and yi for its components. The Born 



Lame coefficients characterize the free energy change of the systems assuming an affine 



displacement of all particles due an imposed external homogeneous linear strain 30|, l31 |. 
Using symmetry considerations it can be readily seen that 

2 

Vb = Xb + -^ (/iB - Pex) , (A7) 

i.e. Ab = {rjB + -Pcx)/2 ~ Vb/'^^ as con&med by the data indicated in Table [B As shown in 
the inset of fig. [151 the Born Lame coefficients per monomer (ffiled triangles) are constant 
for dilute and semidilute densities, but increase strongly for our largest melt densities. 

Stress fluctuations. The two remaining terms Af and fi-p in eq. (IA5|1 characterize the 
fluctuations of the stress tensor 



-Af = pV{^6P,JPyyj (A8) 

-/iF = PV iSP^ydP^y) (A9) 



with Pai3 being the (instantaneous) excess pressure tensor and 5Pq,^ = Pap — {Papj a 
fluctuation. We remind that the mean excess pressure Pex is the averaged trace over the 
instantaneous excess pressure tensor, Pex = ^Tr[PQ,^]^ /d. The above definitions do also hold 
in higher dimensions d > 2 albeit averages may be performed there over equivalent pairs of 



spatial indices. 



related work 



31 



he minus sign in front of Af and /ip is introduced to be consistent with 



33|. By writing the virial W/\^ in the definition of rjp, eq. f|T8|) . in terms of 
the diagonal elements SPaa of the pressure tensor one sees using symmetry considerations 
that 

VF = (-Af) + ^(-/"f). (AlO) 

The reader may verify that eq. flA7|) and eq. flAlOl) are consistent with eq. f|T6|) and eq. flAip . 
Note that eq. (]A4p implies /ib ~ — /ip for low densities where Pex is negligible. While 
Ab = fJ'B by definition for pairwise central forces the Lame coefficients Ap and /ip may differ 
in general. As can be seen from the inset of fig. [T5] (open symbols) or Table [I] we obtain, 
however, to leading order 

Ab/3/p ~ -Ap/3/p ^ -/ip/3/p (All) 
^ (3hll/8 ~ 80 (A12) 

for p < 0.5. In the last step we have used the same reasoning as in eq. f l22|) assuming 
that the bonding potential dominates the moduli for low densities. That eq. (]A11|) holds 



(albeit not rigorously) can again be traced back to the fact that strong two-point (self) 
interaction contributions to rjp (especially due the bonding potential) dominate over the 
distinct correlations ?7F,dist for distinct interactions (/ 7^ /'). Using eq. ( lAlOl) this implies to a 



good approximation — Ap ~ — /ip ~ Vpf^- Coming back to eq. (IA3I) one sees that K is again 
represented by the difference of two large terms Ab and — Ap of same magnitude. Considering 
the statistics of the data visible in the inset of figure, it is clear that our Ap-values do not 
allow to extend the computation of K to smaller densities, just as is was the case using the 
Rowlinson formula. 

Average time- dependent shear modulus. As the hypervirial r/B the Born Lame coefficients 
Ab = f^B reach immediately their asymptotic values (not shown). This is again different for 
the stress fluctuation contribution fip as shown in the main panel of fig. [15] which presents 
the time dependent shear modulus G(t). We use here that according to eq. ( IA2|) and eq. ( lASp 



the shear modulus is given by [32 1 



Git) =fXB- i-Mt)) - Pe.. (A13) 



As it should for a liquid, G{t) vanishes rapidly 33|, |39|. The vertical axis is made di- 
mensionless by tracing the modulus per monomer G{t)P/p. As horizontal axis we use the 
monomer MSD h{t). Interestingly, the data collapse for different densities although we have 
not rescaled the axis with the blob size ^(p) as we did for the compression modulus K{t) 
in fig. [TTl Density effects enter here only through the density dependence of the monomer 
mobility which we have scaled out by using the measured MSD h{t) as coordinate. That 
the scaling of K(t) depends explicitly on the density — and thus according to eq. ( 1A3I) the 
Lame coefficient X^lt) — stems from the fact that the compression modulus couples to the 
density fluctuations which are characterized by ^{p). Not coupling directly to the density 
fluctuations, the shear modulus G(t) and thus f^pit) decay on a local monomer scale. 
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TABLE I: Various properties vs. monomer density p for chains of length N = 1024: the number 
of chains per box M, the root-mean-square bond length /, the root-mean-square chain end-to-end 
distance Re{N) and the gyration radius Rg{N) discussed in sect. IIIlXl the total non-bonded inter- 
action energy e^h per monomer and the interchain monomer interaction energy Cint (Sect. HTlBp . the 
total pressure P (Sect. IIII Cp . the dimensionless compressibility (7t,n = k^Tp/K (Sect. IIIIPp . the 
hypervirial r/B and the excess pressure fluctuation ijp (Sect. HIlEp . and the contributions Ab = PB, 
Ap and pp to the Lame coefficients A and p discussed in the Appendix. The last two columns 
refer to the dimensionless compressibility and the blob size ^ for asymptotically long chains. 
We assume that eq. ^ holds for all densities. Extrapolated values are indicated by stars (★). The 
thermodynamic properties are made dimensionless rescaling them with (3 = l/ksT and the number 
density p. Note that the semidilute blobs become too small around and above p « 0.5 and that 
for = 1024 we have a crossover density p* ^ 0.85/^^/2 0.03, i.e. the semidilute regime ranges 
over about an order of magnitude in density. 
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